Upper Limits on the Continuum Emission from Geminga at 74 

and 326 MHz 



Namir E. Kassim & T. Joseph W. Lazio 

Code 7213, Naval Research Laboratory, Washington, DC, 20375-5351; kassim@rsd.nrl. navy. mil; 

lazio @rsd. nrl. navy, mil 

a^ 
a^ 
a^ 

O : ABSTRACT 

O 

\ We report a search for radio continuum emission from the gamma-ray pulsar Ge- 

^ ■ minga. We have used the VLA to image the location of the optical counterpart of 

Geminga at 74 and 326 MHz. We detect no radio counterpart. We derive upper lim- 
its to the pulse-averaged flux density of Geminga, taking diffractive scintillation into 
O \ account. We find that diffractive scintillation is probably quenched at 74 MHz and 

Q] ■ does not influence our upper limit, S" < 56 mJy (2cr), but that a 95% confidence level 

Q . at 326 MHz is S" < 5 mJy. Owing to uncertainties on the other low-frequency detections 

and the possibility of intrinsic variability or extrinsic variability (refractive interstellar 
scintillation) or both, our non-detections are nominally consistent with these previous 
detections. 
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1. Introduction 

The gamma-ray pulsar Gcminga (PSR J0633+1746, 
IE 0630+178, ICG 195+04) was first identified as a 
compact gamma-ray source in the Galactic anticenter 
(Hermsen et al. 1977; Thompson et al. 1977). Though 
initial speculation included the possibility that this 
source was a pulsar, its identification as such was 
not secure until the detection of X-ray pulsations 
(Halpern & Holt 1992; Bertsch et al. 1992). Since 
pulsars are typically bright in the radio portion of the 
spectrum, intense effort has been devoted to finding 
a radio counterpart to Geminga (Table 1). 

Until recently, only upper limits on the radio emis- 
sion from Geminga could be established. However, a 
number of groups have now reported detecting pulsed 
radio emission from Geminga near 100 MHz (Kuz'min 
& Losovskii 1997a; Kuz'min & Losovskii 1997b; Mal- 
ofeev & Malov 1997; Shitov et al. 1997; Shitov & Pu- 
gachev 1997; Vats et al. 1997; Vats et al. 1999). How- 
ever, subsequent attempts to detect pulsed emission 
at frequencies of 35 and 326 MHz have failed (Ra- 
machandran, Deshpande, & Indrani 1998; McLaugh- 
lin et al. 1999). In an effort to constrain further the 
radio spectrum of Geminga, we have conducted 74 
and 326 MHz observations of it. In §2 we describe 
our observations, in §3 we obtain upper limits for the 
Geminga pulsar that take interstellar scintillation into 
account, and in §4 we discuss our results. 

2. Observations 

Reported previous detections of Geminga have 
been in the range 41-103 MHz (Table 1). Sub- 
sequent attempts to confirm these detections have 
been close to, but outside of, this frequency range. 
The non-detection at 35 MHz (Ramachandran et 
al. 1998) seems problematic, but a combination of a 
low-frequency cutoff in the spectrum and interstellar 
scintillation may explain the detection at 41 MHz and 
the upper limit at 35 MHz. The observations reported 
here are significant in that they are the only attempt 
to confirm the radio emission from Geminga within 
the frequency range of the reported detections. 

Our observations were conducted on 1998 March 9 
with the VLA in the A configuration; Table 2 sum- 
marizes various details. Ionospheric phase variations 
would have vitiated a periodicity search, so our ob- 
servations were designed only to image the pulse- 
averaged emission from the pulsar. The various detec- 
tions suggest that the radio pulse is extremely wide, 



covering perhaps 90° of pulse phase. Thus, our at- 
tempt to detect Geminga should not be affected sig- 
nificantly by searching for a point source at its loca- 
tion. 

Post-processing of low-frequency VLA data uses 
procedures similar to those at higher frequencies, 
though certain details differ.^ At 74 MHz Cygnus A 
served as both the flux density and visibility phase 
calibrator. At 326 MHz the primary flux density cal- 
ibrator was 3C 48 whose flux density was taken to be 
52.2 Jy, and phase calibration was accomplished with 
frequent observations of 4C 14.18. At both frequen- 
cies, after initial phase calibration, several iterations 
of hybrid mapping were used to improve the dynamic 
range. 

Two significant differences for the post-processing 
were the impact of radio frequency interference (RFI) 
and the large fields of view. In order to combat RFI, 
the data were acquired with a much higher spectral 
resolution than used for imaging. Excision of poten- 
tial RFI is performed on a per-baseline basis for each 
10-s visibility spectrum. The large fields of view (11° 
at 74 MHz, 2?5 at 326 MHz) mean that the sky can- 
not be approximated as flat. In order to approach 
the thermal-noise limit, we used a polyhedron algo- 
rithm (Cornwell & Perley 1992) in which the sky is 
approximated by many two-dimensional "facets." 

We did not image the data at the same spectral 
resolution with which they were acquired. Sources in 
our images are identified on the basis of their bright- 
ness. Bandwidth smearing can cause weak sources 
in the outer portions of the primary beam to have a 
brightness below our detection threshold. The side- 
lobes of any such undetected sources contribute to 
the noise in the image. Increasing the spectral res- 
olution has the beneficial result of decreasing the 
amount of bandwidth smearing, however, it also in- 
creases the computational expense and the impact of 
diffractive scintillation. In fact, for reasonable choices 
of the scintillation properties for the line of sight to- 
ward Geminga (§3), the impact of diffractive scin- 
tillation will be unchanged for the range of spectral 
resolutions we have available. The primary factors 
for deciding the spectral resolution are the desire to 
minimize bandwidth smearing vs. computational ex- 
pense. The bandwidths we used — 122 kHz at 74 MHz 



^ A full description of 
data reduction procedures 



low-frequency VLA 
is at (URL: 

http://rsd-www.nrl.navy.mil/7213/lazio/tutorial/). 



Table 1 
Observational Limits on the Radio Counterpart of Geminga Below 1000 MHz 

Frequency Flux Density Reference 

(MHz) (mJy) 

Periodicity Searches 

35 < 100 9 

41 300 7 

61 111 7 

102 < 100 5 
102.5 60 ±95 8 

(5-500) 

102.5 8+^ 6 

103 1000 12 
318 < 1 3 
318 < 3.1 4 
318 < 0.1 11 

325 < 2500 1 
327 < 0.3 9 
327 < 3 10 
430 < 1 3 
430 < 0.6 4 

430 < 0.05 11 

928 < 6 2 

Continuum Searches 

74 < 56 13 

326 < 5.0 13 



References.— (1) Mandolesi et al. (1978); (2) Seiradakis (1981); (3) D'Amico (1983); (4) Fauci et al. (1984); 
(5) Kuz'min & Losovskii (1997b); (6) Shitov & Pugachev (1997); (7) Shitov et al. (1997); (8) Malofeev & 
Malov (1997); (9) Ramachandran et al. (1998); (10) McLaughlin et al. (1999); (11) Buderi et al. (1999); (12) Vats 
et al. (1999); (13) this work 

Note. — Only the work of McLaughlin et al. (1999) and that reported here have taken interstellar scintillation 
into account explicitly in setting upper limits. 









Table 2 
VLA Observing Log 






V 

(MHz) 


B 

(MHz) 


Bi 

(kHz) 


T Beam 

(hr) (" X ") 


A/ 
(mjy bcam^^) 


F 

(mJy) 


74 
326 


1.5 
3 


122 
684 


5.4 25 X 23 @ 64° 
5.4 5.1 X 4.8 @ 75° 


28 
1.4 


-13 
0.28 



Note. — B is the total observing bandwidth while Bj is the bandwidth used in imaging. T is the on-source 
integration time. A/ is the rms noise level in the image, and F is the measured flux density at the location of 
Geminga. 



and 684 kHz at 326 MHz — gave modest spectral reso- 
lutions (a few), modest bandwidth smearing near the 
edge of the primary beam (~ 10%), and reasonable 
imaging times (less than a few days of computation). 

Figure 1 shows our 74 MHz image of the immedi- 
ate region around Geminga. The off-source rms noise 
level on this image is 28 mJy beam"-'^, and the es- 
timated thermal limit is 25 mJy beam~^. Figure 2 
shows our 326 MHz image of the immediate region 
around Geminga. The off-source rms noise level on 
this image is 1.4 mJy beam~^, and the estimated 
thermal limit is 0.5 mJy beam~^. 

One additional impact on the 74 MHz observa- 
tions is ionospheric refraction, which shifts the field 
of view without distorting the brightness distribution 
(Erickson 1984). Shifts at 74 MHz are typically a 
few arcminutes and vary on time scales of tens of 
minutes. Self-calibration "freezes" out this refraction 
but destroys absolute position information. Fortu- 
nately, the primary beam is large enough that a typ- 
ical image contains many tens of sources identified in 
the 1400 MHz NRAO VLA Sky Survey (Condon et 
al. 1998). We constructed a reference grid using 42 
NVSS sources that were either unresolved or slightly 
resolved (diameters < 30"). In determining the ref- 
erence grid, we checked for any systematic trends as 
a function of position within the 74 MHz image, flux 
densities of the sources, or angular diameters of the 
sources. None were found. We shifted the 74 MHz 
image by the median offsets, 3'.'9 in right ascension 
and — 58'.'4 in declination. Based on the variances of 
the offsets from the reference sources, our 74 MHz 
astrometry should be accurate to approximately 5", 



a fraction of the beam diameter, for an individual 
source. Figure 1 has been corrected for refraction. 

3. Upper Limits and Interstellar Scintillation 

In this section we address our upper limits on Ge- 
minga's flux density and the impact of scintillation 
on them. We shall follow the procedure described by 
McLaughlin et al. (1999) for dealing with the effects 
of interstellar scintillation (ISS). 

Briefly, their procedure involves convolving the 
known probability density functions (pdfs) for the 
diffractive ISS gain and noise in the image to produce 
a pdf for the measured flux density for the scintil- 
lating source. This pdf is then multiplied by a prior 
for the intrinsic flux density of the source (Geminga) 
to form the posterior for the intrinsic fiux density of 
the source. As our images were constructed using an 
FFT, without specifying a zero-spacing fiux density, 
we have assumed (and verified) the noise in the image 
has a zero-mean Gaussian pdf. 

In addition to the rms noise in the image (Ta- 
ble 2), we require two additional parameters: a mea- 
sure of the strength of diffractive ISS and a mea- 
sured fiux density at the location of Geminga. The 
former is provided by TVdisSi the number of "scin- 
tilles" or scintillation intensity maxima occurring dur- 
ing the entire observation (Cordes 1986; McLaugh- 
lin et al. 1999). For Geminga, there are no pub- 
lished measurements of the relevant ISS parameters, 
the scintillation time At^ or the scintillation band- 
width Avd- McLaughlin et al. (1999) estimated these 
parameters to be At^ ~ 275 s (z^/327MHz)2-2 and 



Ai/rf « 1.5MHz(j//327MHz)4-4. Their values im- 
ply a scattering measure toward Gcminga of SM ~ 
10"'*'* kpc m~^°/^. While consistent with the local 
distribution of scattering material, there are also large 
variations in the value of SM for pulsars at distances 
comparable to that of Geminga: PSR J0437— 4715 
with D = 0.18 kpc and SM = lO^^.s y.^^ ^1-20/3 
(Sandhu ct al. 1997; Johnston, Nicastro, & Koribal- 
ski 1998) vs. PSR B0950+08 with D = 0.125 kpc 
and SM = lO-''-^ kpc m-^o/s (Qwinn ct al. 1986; 
Phillips & Clegg 1992). We shall adopt McLaugh- 
lin et al.'s (1999) value for SM but note how a dif- 
ferent value affects our results. For our integration 
time and receiver bandwidth, the resulting values are 
A^Diss == 910 at 74 MHz and A^diss = 15 at 326 MHz. 
We obtain the flux density at the location of Ge- 
minga by integrating over an area of 1 beam at the 
position of Gcminga. We use the term "flux density 
at the location of Geminga" because of the noise in 
the image. As the noise is zero-mean, there are both 
positive and negative intensity variations. A weak 
source combined with a large negative noise deviation 
could result in a measured flux density quite different 
than that of the intrinsic flux density of the source. 
Given the sizes of the beams (6" at 326 MHz and 25" 
at 74 MHz), this integration also takes into account 
the uncertainty in the position of Geminga which is 
less than 1" (Bignami et al. 1993; Caraveoet al. 1996). 

3.1. 74 MHz 

We begin by determining the upper limit on Ge- 
minga's flux density if ISS was not an issue. The po- 
sition of Geminga's optical counterpart, even taking 
into account its proper motion, is known to accuracy 
far better than the beam diameter of 25". We there- 
fore adopt an upper limit of 56 niJy, twice the rms 
noise level. 

As Ndiss -^ 00, DISS becomes more and more 
heavily quenched, and the pdf of the DISS gain ap- 
proaches a (5-function centered on unity. We con- 
sider A^Diss ^ 10"^ sufliciently large that DISS will 
not affect our upper limit. If our adopted value 
of SM were too large, and the correct value were 
SM ^ lO^'''"' kpc m^^^/-^ — comparable to that for 
PSR B0950-h08— then TVdiss = 44. We would ob- 
tain an upper limit for Geminga's flux density of 
S' < 86 mJy at the 95% confidence level. 



3.2. 326 MHz 

If we ignore ISS, the 6" beam remains large enough 
that wc adopt an upper limit of 2.8 mJy, again twice 
the rms noise level. 

Our adopted value of SM implies iVniss — 15. 
In contrast to the situation at 74 MHz, DISS can 
have a substantial impact on the adopted flux density 
at 326 MHz. We find an upper limit on Geminga's 
flux density of S" < 5.0 mJy at the 95% confidence 
level. If our adopted value of SM were too large, and 
the correct value were SM ^ lO^'"'"^ kpc m^^^/"^, then 
A^Diss = 5. We would obtain an upper limit for Ge- 
minga's flux density of S* < 13.5 mJy at the 95% 
confidence level. 

4. Conclusions 

Figure 3 incorporates our upper limits with the 
limits and detections of Geminga below 1000 MHz in 
the literature; Tabic 1 summarizes the data shown in 
Figure 3. 

Can we reconcile our nondetection with the detec- 
tions near 100 MHz? One simple possibility is that 
the pulse-averaged radio emission from Geminga is 
intrinsically time variable. In this case, our nonde- 
tection is quite easy to reconcile with the various de- 
tections. Vats et al. (1999) argue that Geminga is 
intrinsically time variable in order to explain their 
detection with a fiux density of 1 Jy at 103 MHz, and 
they speculate on intrinsic mechanisms that could be 
responsible. 

Even if we treat the Vats et al. (1999) detec- 
tion at 103 MHz as an outlier, though, our nonde- 
tection is just consistent with the other detections 
near 100 MHz. Shitov & Pugachev (1997) report a 
flux density of 8+2 mJy at 102.5 MHz. Clearly, if this 
is the intrinsic flux density, then there is no difficulty 
reconciling this with our nondetection. 

There are two difficulties with assuming a flux den- 
sity of 8 mJy at 102 MHz, though. First, it requires 
Geminga to have an extremely steep spectrum be- 
tween 61 and 102 MHz, a « —5, though, Geminga's 
radio emission mechanism might be different than 
that of other pulsars. Second, it is difficult to explain 
other detections at 102 MHz. Some combination of 
intrinsic and ISS-induced intensity variation could in 
principle explain this dynamic range, but, in practice 
it might require improbably large DISS gains. 

A more likely value for the 100 MHz flux density is 



near 50 mJy, in which case Shitov & Pugachev (1997) 
detected it in an ISS-depressed state. This explana- 
tion is appeahng in two respects. First, a DISS gain 
of 0.1 would depress the flux density from near 50 mJy 
to near 5 mJy, and a DISS gain this small has a prob- 
ability of occurrence near 90%. Second, the spectral 
index between 61 and 102 MHz would be —1.6, more 
characteristic of radio pulsars. In this case, the ex- 
pected flux density at 74 MHz is 80 mJy, above our 
upper limit for our nominal value of ^VdisSj but below 
it if our adopted value of SM is too large. 

In reconciling our nondetection at 74 MHz with the 
detections, we have relied on the robustness of the 41 
and 61 MHz observations. In reporting their detec- 
tions at these frequencies, Shitov et al. (1997) did 
not indicate to what extent DISS had been quenched. 
As DISS is more likely to depress rather than enhance 
the flux density of a source, the 41 and 61 MHz detec- 
tions could represent lower limits to Geminga's flux 
density. If so, reconciliation would be more difficult. 

One aspect of ISS not addressed quantitatively is 
the impact of refractive ISS. At low frequencies, ISS 
has an ancipital nature: DISS occurring on "short" 
time scales and RISS occurring on "long" time scales. 
RISS is more difficult to quench than DISS because its 
time scales are typically hours to days, so the flux den- 
sity of Geminga could be enhanced or suppressed for 
an entire observing run. RISS could serve as the ex- 
trinsic mechanism that would contribute both to our 
non-detection (and that of Ramachandran et al. 1998 
at 35 MHz) and the detections near 100 MHz. 

RISS is more difficult to handle quantitatively be- 
cause its pdf is not as well constrained as the pdf for 
DISS. Some indication of the effect of RISS can be ob- 
tained by adopting a log-normal pdf (Ishimaru 1978, 
Chap. 20). Gupta, Rickett, & Coles (1993) found 
a refractive intensity modulation variance m^ « 0.2 
for PSR B0950-h08 at 74 MHz. This variance, when 
used with a log-normal pdf, implies a (95% probabil- 
ity) range of RISS-induced variations of a factor of 
a few. Thus, these low-frequency measurements and 
upper limits could have an additional factor of two or 
so uncertainty. 

In summary, we have not detected Geminga at 
74 MHz {S < 56 mJy) or 326 MHz {S < 5 mJy). 
We can reconcile our nondetection with previous low- 
frequency detections by invoking intrinsic variability 
or extrinsic variability (refractive interstellar scintil- 
lation) or both. 
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Fig. 1. — The Geminga field at 74 MHz. The cross marks the location of the optical counterpart of Geminga 
(Caraveo et al. 1998); the uncertainty in the location of Geminga is approximately 100 times smaller than the size 
of the cross. The beam is 24'.'6 x 23'.'3 at a position angle of 64°, and the rms noise level is 28 mJy beam~^. Contour 
levels are 28 mJy beam^^ x —2, 2, 3, 5. Also indicated are two sources from the NRAO VLA Sky Survey (Condon 
et al. 1998). 



Fig. 2. — The Geminga field at 326 MHz. A cross marks the location of the optical counterpart of Geminga 
(Caraveo et al. 1998); the uncertainty in the location of Geminga is approximately 100 times smaller than the size 
of the cross. The beam is 5" x 4'.'8 at a position angle of 75°, and the rms noise level is 1.4 mJy beam~^. Contour 
levels are 1.4 mJy beam^^ x —3, 3, 3, 5, 7.07, 10, 14.1, .... Also indicated are a number of sources from the NRAO 
VLA Sky Survey (Condon et al. 1998). 



Fig. 3. — The radio spectrum of Geminga below 1000 MHz. The upper limits reported in this paper are marked 
with a star. Limits and detections are summarized in Table 1. With the exception of the upper limit by McLaughlin 
et al. (1999; < 3 mJy at 326 MHz), all other upper limits have not considered the effects of interstellar scintillation. 
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